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Abstract 

A phase space coalescence description based on the Wigner— 
function method for cluster formation in relativistic nucleus-nucleus 
Q ' collisions is presented. The momentum distributions of nuclear clus- 

^ . ters d,t and He are predicted for central Au(11.6AGeV)Au and 

Si(14.6AGeV)Si reactions in the framework of the RQMD trans- 
port approach. Transverse expansion leads to a strong shoulder- 
arm shape and different inverse slope parameters in the transverse 
spectra of nuclear clusters deviating markedly from thermal dis- 
tributions. A clear "bounce-off" event shape is seen: the aver- 
aged transverse flow velocities in the reaction plane are for clusters 
larger than for protons. The cluster yields -particularly at low p t 
at midrapidities- and the in-plane (anti)flow of clusters and pions 
change if suitably strong baryon potential interactions are included. 
This allows to study the transient pressure at high density via the 
event shape analysis of nucleons, nucleon clusters and other hadrons. 

1 Introduction 

One of the challenges of modern heavy ion physics is the extraction of the 
equation of state and transport properties of extremely dense and excited nu- 
clear matter. In particular, the study of matter at high net baryon density 
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has received much attention recently. QCD - as the accepted theory of strong 
interaction - contains chiral symmetry (in the limit of massless quarks) which 
is spontaneously broken in its ground state, the QCD vacuum (see e.g. lattice 
calculations H ) . A rapid restoration of this symmetry with increasing baryon 
density is predicted by all approaches which embody this fundamental aspect 
of QCD ||, ||. Nucleus-nucleus collisions in the bombarding energy region of 
baryon stopping may therefore be favourable to study such medium effects as 
compared to ultrahigh energies, for which the two colliding nuclei may become 
transparent. Beam energies between 10 to 15 AGeV - as studied experimen- 
tally at the BNL-AGS d]-f§] ~ seem to be well suited to stop two heavy ingoing 
nuclei and to create the desired high baryon densities. This has been shown 
by transport calculations based on hadronic excitations and rescattering like 
the RQMD approach (strings, resonances) §, 0, the ARC- @] or the ART- 
model [[[3) (resonances). The observation of stopping in the AGS experiments 
has been unclear for quite some time. However, all experimental groups now 
confirm [Q-(t]] the predicted large baryon stopping in central collisions || Q. 

An observable consequence of the formation of dense nuclear matter - far be- 
yond the groundstate - is the emergence of collective flow driven by compression- 
induced pressure jlj|-[jl8|. Mean fields |l5| may give important contributions 
to this pressure and could therefore be accessible to experimental observation, 
just as in the lGeV region [|l9|. The bounce-off for protons has been observed 
at 10 GeV/n pQ] as well as azimuthally asymmetric particle correlations in the 



projectile hemisphere 21 . These experimental discoveries encourage us to in- 
vestigate the formation of nuclear clusters - as compared to light hadrons - for 
which flow can even dominate the momentum spectra [ p~7| . 



We follow our earlier work on deuterons 22 and extend the phase space 



coalescence picture to light clusters with A < 4. The basic ingredients of the 
cluster coalescence, i.e. the source function provided by the RQMD and the 
parametrization of cluster wave functions are described in Sect. 2 and 3. The 
parametrization of baryonic mean fields is described in Sect. 4. In Sect. 5 par- 
ticular features in the momentum distributions of nuclear clusters are discussed 
as signatures for flow and event shape correlations at the hadronic freeze-out. 
In order to demonstrate sensitivities to baryonic mean fields final observables - 
rapidity distributions, mt-spectra, directed Row p x (y)- are compared for two ex- 
treme scenarios: one with a density-dependent quasi-potential between baryons, 
and the other without (cascade). 

It should be mentioned that results for cluster yields have been calculated 
based on the thermal model and on the single-nucleon momentum distributions 
[H |J, II H H> HI HI HI H H HI- These results predict that the spectra 



of clusters and nucleons have essentially the same shape. Here, we are going 
to demonstrate that flow invalidates the basic assumptions underlying these 
simple models. In turn, we can use the amount of 'scaling violation' of cluster 
spectra as compared to proton spectra to assess the strength of collective flow 
in nucleus- nucleus reactions. 
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2 Coalescence of Clusters in Phase Space 



We combine a dynamical description of the first violent stage in nucleus-nucleus 
reactions from the RQMD model with a cluster formation model which is based 
on the single-particle phase space distributions at freeze-out. RQMD is a semi- 
classical transport theoretical approach and does not take into account forma- 
tion of nuclear bound states (e.g. deuterons) dynamically. However, the small 
binding energies and the associated quantum mechanical formation time from 
the uncertainty principle suggest that nuclear clusters are produced mainly much 
later after the violent interactions have ceased (freeze-out), i.e. cluster formation 
rates can be calculated from the nucleon distributions at freeze-out. In order 
to calculate light nuclear cluster distributions (for A<4) we use the Wigner- 
function technique in phase space. This phase space coalescence approach was 
already applied to deuteron production at bombarding energies around lAGeV 
@ ||, 10-15AGeV H |6[ and 160-200AGeV §|. 

The validity of the combined RQMD+coalescence approach for the cluster 
formation clearly depends on that the transport model describes the dynamical 
evolution up to the freeze-out stage reasonably well. The relativistic quantum 
molecular dynamics approach (RQMD 1.07) |^| employed for the calculations 
presented here combines the classical propagation of particles with excitation of 
hadrons into resonances and strings. Secondaries (emerging from the decaying 
resonances and strings) undergo subsequent interactions, both with each other 
and with the ingoing baryons. Note that the RQMD-results compare well with 
experimental single particle and 2-body correlation data || [?], [h], |llj . 

In the following, we use the non-relativistic Wigner-function formalism which 
may be justified for nuclear clusters in view of the small binding energies. The 
formation rates are calculated at equal time in the common rest frame of the 
corresponding cluster nuclcons immediately after their freeze-out. Having this 
in mind, we suppress the explicit time dependence and reference to the chosen 
Lorentz frame in all following expressions. Note, however, that all results pre- 
sented in this work include an implicit integration over all freeze-out times and 
Lorentz transformations back into the original observer system. 

The Wigner-function for a single particle 

J d 3 y < x + y/2\V >< V\x-y/2 > e~^ M (1) 

is the closest analogue to a probabilistic distribution in phase space which one 
can get from quantum mechanics. Therefore, its identification with the phase 
space distribution fa as p w has been frequently employed in semi-classical 
calculations. Neglecting the hopefully small effects from binding energies, the 
formation probability for a cluster can be expressed as an overlap integral be- 
tween the Wigner-function which corresponds to the cluster wave function and 
the N-body nucleon phase space distribution at freeze-out |?9|, |34|, |3£ , 
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The N-body phase space distribution has to be constructed from the single par- 
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ticle "source function" which is defined by the "freeze-out" positions x^ 
momenta of nucleons after their last scattering or decay. 
The Wigner-density of an M-nucleon state has the form 



and 



[PMl 



E 

[TT 3 ]x[SS 3 ] 



(\SS 3 >< SS 3 



\TT 3 < TT 3 \) 



(2) 



x W[tt 3 ] x [ss s ] {xx,pi;...,x M , Pm) 

with the normalization tr^Af] = 1. The product [TT 3 ] x [SS 3 ] denotes the 
set of all 2 2M possible internal couplings with proper total isospin T, T 3 and 
spin S,S 3 . Note that the phase space part W\tt 3 ]x[SS 3 ] has to provide the 
correct symmetry concerning particle exchange to ensure that all states are 
totally antisymmetric. In the semiclassical approximation it is assumed that 
the Wigner-function does not contain dynamical correlations with respect to 
spin and isospin. We therefore employ the statistical assumption and assign all 
many-nucleon states which are allowed from the Pauli principle the same weight 
for a given position and momentum distribution: 

We consider only spin averaged Wigner-densities W^[tt 3 ] x [SS 3 ] ~ 1 /2 M Wpyr 3 ] ■ 
Furthermore, the coupling of M particles to a given total isospin T is assumed 
to be equal 

W [T t 3 ] x [ss 3 ] » ITm Wt s = 9 Wt 3 (3) 



( M z) 2M 



In this approach Wt 3 contains all spin states and all ( z ) states in different 
isospin-multiplets for a M-particle combination (M, Z) with given charge Z = 
T 3 + M/2. 

In order to approximate Wt 3 we use the RQMD model that provides the 
phase space distribution of nucleons with given isospin. We identify Wt 3 with 
the product of single particle distributions 
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t 3 



N M 



where N = iV„ 



_i=Z+l 

N n is the total number of nucleons and 



M 



J] {2TTh) 3 f n (x i ,p i ) 



(4) 



N p := / d 3 xd 3 p f p (x,p) , N n := / d 3 xd 3 p f n (x,p) 



(5) 

Eq. U can be interpreted as a statistically uncorrelated emission. It defines the 
probability-density to find a given nucleon combination (M,Z) in certain phase 
space regions. Inserted in Eq. ^3] it fulfills by construction the trace normalization 
tr[pAf] = E 

The cluster wave function which is assumed to be non-relativistically here 
factorizes into a collective and a relative part 

\*c(P) >= 7 ^ WH ei p -*<f>(t 1 ,...,t M -i)\SS 3 > \TT 3 > (6) 



(27rS) 3 / 2 



4 



where X = (xi + ...+xm)/M and P = pi + ...+pM. The U(xi, xm) (i=l,...,M- 
1) are the M — 1 relative coordinates of the relative cluster wave function (p. 
S, S3, T, T 3 are the spin and isospin quantum numbers of the cluster state. The 
Wigner-density of the wave function in relative coordinates is defined by the 
Wigner-transformed projection operator 

p%(t u ft; ? u -i,fa-i) ■= \TT 3 > \SS 3 > p% < SS 3 \ < TT 3 \ (7) 

with 

9c : = J 01 + -,tM-i + yM-i/2)<f>*(.h - j/i/2, -Jm-i - Vm-i/2) 
Xe - *** ... e-i^-iw-i^ ... d 3 y M -i (8) 

The formation of cluster states is finally determined by the trace over the source 
density pm and the projector on the individual cluster wave function >< 



[pm] W (xi,pi; .■■\&M,fM)Pcftfi><lii -;*m-i>3m-i) 



X6 3 {P- (pi + ...+PA/)) 



dx\d 3 px d 3 XMdp 



(2tt?i) 3 "' {2ith)- 



M (9) 



The absolute number of states is obtained by multiplying Eq. || with the total 
number of M-nucleon states and summation over all possible spin states 
N$ ■ Inserting Eq. || and || the semi-classical coalescence formula reads finally 



dN 

d^ = 9NS 



N\ (M 
M 
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dx\dpl...dxl I dpl I 8 A {P - {pi 
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Pm)) 



J| fn(Xi,Pi) 



U=l J Li=Z+l 

3u D n-^A 4t 



9c (h,qi;-;t M -i,qM-i) (10) 



For d, t, ''He and He states the momentum distributions are explicitly given 

by 



dN(d) 
d 3 P 



/(d) N s (d) 
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N 2 



dx 1 dp 1 dx 2 dp 2 f n (xi,pi)f p {x 2 ,P2) 



(11) 
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X ~j^3 f dx l d Pl--- dx 3 d Plfn(xi,Pl)fn{x2,P2)fp(x3 : P3) 

Xpf 5 (P (12) 



d3p g i*He)N a ?Be) 



X Jp f dx l d Pi--- dx l d Plfn(xi,Pl)fp(x2,P2)fp(x3,P3) 

x Pm c 6 (P - (Pi +P2+P3)) (13) 



<PP av ' ov J \ A) \2 

1 /" 

X Jri dx l d pl--- dx 4 d Plfn( x l,Pl)fn{x2,P2)fp{x3:P3)fp( x i,P4) 

x p 4 ^ c 5 (p - (pi + p 2 + p 3 + p 4 )) (14) 

In a Monte Carlo formulation - appropriate for the application to microscopic 
transport calculations - these formation rates can be expressed by the general 
coalescence formula for M-body cluster 

dN M = gN s (^ ^2 Pc (4.fe;-;*iM-i)9iM-i))\ 

n<~<*M 

xd 3 t ll d 3 q ll ...d 3 U M _ 1 d 3 q lM _ 1 . (15) 

< ... > denotes event averaging. The sum runs for each event over all M-nucleon 
combinations. Note the necessary condition i\ < ... < im which prevents the 
double counting of equal particle pairs. The coordinates in position and mo- 
mentum space are taken at equal time in the M-nucleon rest frame (i.e. P = 0) 
immediately after all cluster nucleons have frozen out. The calculated num- 
bers contain higher mass fragments by construction. The number of A> 4 
clusters is small, however, for rapidity values \y — y m id\ < 1- The factor g 
contains spin and isospin projection as described above. After the summation 
over all possible spin states the statistical corections are g(d)Ns(d) = 3/8, 
g(t)N s (t) = g( 3 Re)N s ( 3 Re) = 1/12 and 3 ( 4 He)N s ( 4 He) = i/g 6 _ Feed down 
effects from the production of excited t and He states are expected to be small 
(< 15%, Q) and will be neglected in the present studies. 

The statistical approximation employed here is expected to break down in 
regions where the binding energies and the quantum dynamics play an essential 
role, e.g. in case of spectator matter fragmentation. Deviations from the statis- 
tical limit could give further insight into the fermionic (quantum) dynamics of 
the many-body system and final state effects like e.g. Coulomb distortion. 
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3 The Parametrization of Cluster Wave Functions 



For the deuteron we assume a Hulthen-wave function derived from a Yukawa- 
type potential interaction |^2|. Ea] 

x ,-, _ | [exp(-a|xi - x 2 |) - exp(~6|fi - x 2 \)]- (16) 
Fi - x 2 \ 

In order to get a simple expression for the Hulthen Wigner-density the wave 
function is approximated by a sum over 15 centrally symmetric Gaussian wave 
packets 



*(r) = > >G(r) = Va, ( ^ ) exp(-c,r 2 ). (17) 



2f; 



3/4 



The Wigner-density of this sum can be calculated analytically 
PD{r,q) = 8^<z?exp(-2cir 2 )exp ( — 



+ 16 a,aj 



4c,-c 



i 1 j 



3 V(c« + Cj ) 2 



3/4 



2 \ ( ~q 2 \ ( n a 



x exp I — r I exp I I cos I 2 -rqcos9 rq 1(18) 

\a + Cj J \ a + Cj J \ a + cj ) 

with q = \px — P2I/2 , r — \x\ — x-z\ as the relative position and momentum 
coordinates. 

For triton, 3 He and 4 He states we use 3- and 4-particle harmonic oscillators 
with different coupling strength. Such an approximation has been used already 
in momentum coalescence studies [^4], |27], and is well known in nuclear 

structure physics (see e.g. @). We adopt such a form of the wave function 
because of the separability in collective and relative motion even on the level of 
M-particle states. Moreover, a harmonic M-particle wave function can always 
be written as a product of single particle oscillators which leads to a Wigncr 
density of the form 

M-l 

q = S(q M - (qi + ... + q M -x)) I] 8e-ft ^ (19) 

j=l 

The tj — ^2(C)jiXi are given by the linear transformation of the original carte- 

i 

sian coordinates C . The generalized relative momenta are q 3 ■ = (C _1 ' + ) k Pk- 

k J ' 
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For complicated systems the transformation C can contain particle masses and 
different coupling constants. Our purpose has been a simple parametrization 
rather than taking into account detailed information in the wave function as 
e.g. Coulomb repulsion. Therefore, we use only one coupling constant D for 
each cluster and an equal mass of protons and neutrons m which leads to a 
representation in so-called Jacobi coordinates: For t and 3 He states the two 
relative coordinates are t\ = ^3D/2(x 2 — x\), t 2 — V2D(x 3 — (x\ + X2)/2) 
and Xq = Xq = (3D/2m) 1 / 4 while for 4 He states t\ = \/2D{x 2 — X\), t 2 = 
y/8D/3{x 3 - (x 1 +x 2 )/2), t 3 = V3D{xi - {xi + x 2 +x 3 )/3) and = x\ = X% = 
(2D/m) 1 / 4 . The coupling constants are adjusted to the mean square charge 
radii of the diverse cluster states (see Table 1 and [f^] ) 



4 Baryonic Mean Fields 



At nuclear ground state density the nuclear mean field may be decomposed into 
two large pieces: an attractive scalar field provided by the quark condensate 
and/or correlated two-pion exchange (the a field), and a repulsive vector poten- 
tial (the u> field) |Q which is in accord with Dirac phenomenology for optical 
potential calculations Q-Q in p+A studies and QCD sum-rule estimates 
Not much is known about the strength of the mean fields at large net baryon 
densities and temperatures predicted in all present transport approachs for the 
ultrarelativistic regime. It is expected that the momentum dependence |35, [5l| , 
the excitation into resonances [55|] and the transition to quark matter [56| will 
play a crucial role for the created mean fields. Several new ideas are currently 
under active investigation: Medium properties of hadrons (e.g. of the u) meson 
which is responsible for vector repulsion |^]) or quark and gluon condensates, 
which break the approximate scale and chiral symmetry of the QCD Lagrangian 
in the vacuum, could modify the scalar field essentially |E6l || . 

In the following we demonstrate the sensitivity of flow observables to mean 
fields by comparing two schematic cases: In the first case the potentials are 
switched off (i.e. the cascade mode is used). The second scenario uses potential 
type interactions which define effective baryon masses in a medium Q 



Pi 



Vi = 



and thus simulate the effect of mean fields. Here 

Pij\ , P 



(2m N )' 1 V l =+lYl 



7+1 



E 



8li 
Pa 



(20) 



(21) 



with pij a Gaussian of the CMS distance vector normalized to one, po ground 
state matter density and a = — 0.4356GeV, (3 = 0.385GcV, 7 = 7/6 parameters 
which are adjusted to the saturation properties of nuclear matter (binding en- 
ergy and compressibility) . The parameter fit was done by taking the expectation 
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value of the total energy per nucleon for ideal gas (plane) wave- functions taking 
into account antisymmetrization. A Hamiltonian can be constructed which con- 
serves the effective mass shell constraints of Eq. 20 [pi and is employed for the 



dynamical evolution in RQMD. For the presented calculations we have chosen 
to let the particles interact at equal times in the global equal-speed system of 
projectile and target. 

As has been stated in |57|, the experimental data for nucleus-nucleus reac- 
tions at 10-15AGeV indicate more repulsion than just given by a pure density 
dependence as in Eq. (3). In fact, it was found that the quasi-potentials do not 
affect the final distributions at all, if their strength is below a 'critical threshold'. 
This additional repulsion is probably caused by the momentum dependence of 
nuclear forces. In order to explore the possible role of mean fields, in Ref. |37j 
we hardened the density dependence of the potentials until we achieved agree- 
ment with experimental proton spectra which were available at that time. The 
attractive 2-body force in the AA and NB* channel have been switched off - 
qaa = «nb* = - (thus explaining the index pair (ij) in Eq. (3)). Here, we 
use the same potentials. We note that this approach has predictive power for 
the cluster spectra, because the strength of the potentials and thus the amount 
of flow has been fixed from the proton spectra alone. 

What does the existence of a 'critical threshold' above of which only the 
collective forces 'win' against the randomization of the motion by stochastic 
collisions and decays mean? It indicates that the mean field effect cannot be 
isolated from the other - stochastic - interactions which are present in the 
system. For instance, the initial baryon density is essentially fixed solely from 
the stochastic interactions, because the degree of energy degradation which the 
ingoing nucleons experience is mostly determined by the multiple collisions with 
nucleons from the other ingoing nucleus and with secondaries. More recent 
versions of RQMD contain somewhat more realistic interactions than used here 
for the presented calculations (see Ref. j5£| for a discussion). For instance, the 
assumption of istropic heavy baryon resonance decay leads to an overestimation 
of the nuclear stopping power as it was noted recently in Ref. [p8|| . A detailed 
study of the interplay between the effect of collisions and mean fields based on 
RQMD is currently undertaken by one of us J5?|. On the other side, we do 
not expect that our qualitative conclusions about the effect of mean fields on 
cluster flow will be reversed by more realistic calculations. In fact, a smaller 
initial baryon density means that the mean fields have to get stronger to achieve 
the same amount of collective flow. 



5 Results and Discussion 

The production of clusters has recently been measured and analyzed for central 
and peripheral reactions p+A, Si+A and Au+Au at AGS-energies |30| [H], ^, |62[ 
and S+A at 200AGeV J63[ . Comparisons between coalescence calculations for 



deuterons at pt — and measurements for pA reactions have also been discussed 
|. We will first show that the phase space coalescence in combination 
with the RQMD freeze-out describes absolute values and momentum distribu- 
tions of deuterons in accordance with measurements for various nucleus-nucleus 
collisions. In Fig. [l]we compare our recently published results for deuterons in 
central Si+A reactions at 14.6AGeV |§| with E802 data 0. In Fig. | cal- 
culations for transverse mass spectra of protons and deuterons in the reaction 
Au(11.6AGeV)Au (b<3fm) at rapidity y — 1.3 are compared to preliminary 
E866-data j(34|] for central events. A comparison between calculations and data 
requires proper event selection according to experimental trigger conditions and 
acceptance corrections for the theoretical calculations (65). On the level of the 
systematic errors in the measurements (« 15%, ||| ^l|), however, we find good 
agreement - even for the strong slope parameter splitting between protons and 
deuterons in massive reactions. 



5.1 Transverse Expansion and Cluster Flow 

The formation of transverse blast waves was first proposed in Ref. |66| where 
pion and proton transverse momentum spectra around lAGeV incident beam 
energy were analyzed. The most prominent observables are the characteristic 
shoulder-arm shape and different apparent temperatures for particles with dif- 
ferent mass. These are caused by the overlay of rather small local momenta and 
collective motion which have been produced during the expansion phase of the 
reaction. Several investigations followed in the low energy regim e O, |67 , |(38f 
as well as for ultrarelativistic nucleus- nucleus reactions |6^, |7^, |7l|, |72|, |24 3f|. 

In the following we will discuss the momentum spectra of nuclear clusters 
which show a strong dependence on such phase space correlations. In partic- 
ular heavy clusters like 4 He can serve as a very promising tool to determine 
the phase space picture as e.g. provided by the microscopic transport calcula- 
tions. The strongest flow and mean field effects are achieved in central reac- 
tions Au(11.6AGeV)Au. Variations with the reaction size and life time of the 
high density zone are studied by comparison with results for Si+Si reactions at 
14.6AGeV. All results presented here have been calculated for central impact 
parameters (Au+Au: b<3fm, Si+Si: b<lfm). The weak decay of hadrons after 
the freeze-out has been suppressed. 

Rapidity distributions and transverse momentum spectra of p, d, t and 4 He 
cluster are shown in Figs. || and [| for central Au+Au (11.6AGeV) and Si+Si 
(14.6AGeV) reactions. The figures contain cascade (solid histograms) and po- 
tential calculations (bold solid histograms). The solid lines in Fig. ^ show 
Boltzmann distributions with temperature parameters adjusted to fit the trans- 
verse spectra for p t >2GeV/c (Au+Au) and p t >lGeV/c (Si+Si) as calculated 
with baryon potential-interaction. 

In Au+Au collisions all rapidity distributions peak at midrapidity indicating 
strong stopping in accordance with earlier predictions and preliminary data for 
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protons [22| . The transverse distributions have a strong shoulder-arm shape 
which deviates markedly from distributions expected from a purely thermalized 
fireball. The shoulder-arm shape becomes most prominent for heavy clusters. 
For He clusters a peak even appears at finite pt- The high momentum tail of 
the transverse spectra exhibit different 'apparent' temperatures for clusters with 
different mass while a thermal system would predict similar slope parameters 
[j§H . Note that the extracted temperature values depend strongly on the p t -cut 
choosen. The absolute values extracted by exponential parametrizations always 
lead to additional systematic errors in the absolute yields according to our cal- 
culation overestimating the cluster yields substantially. Therefore, comparisons 
of rapidity distributions between calculations and extrapolated data must be 
done very carefully. 

In contrast, all rapidity distributions are essentially flat for the light system 
Si+Si. The transverse spectra are in good agreement with a thermal fit and 
show temperature parameters which are almost equal for all states. 

The characteristic deviations from thermal distributions are caused by strong 
transverse expansion and collective flow particularly in massive reactions like 
Au+Au. The flow-correlations at the microscopic freeze-out are shown in Fig. 
U which contains calculations for the freeze-out velocities and density profiles 
of protons, deuterons, tritons and 4 He. The velocity profiles for all clusters are 
similar. They exhibit a convex shape and saturate at «0.7c for both reactions. 
In the massive system Au+Au the freeze-out densities have a complicated shape 
which peakes around 5fm. Most of the nucleons freeze out at larger distance. 
This is indicated by an average freeze-out radius of w lOfm (see Table 2). The 
strong transverse expansion in Au+Au collisions is caused by the considerable 
baryon stopping and the pile up of high particle densities near to the reac- 
tion center. During the expansion phase comoving particles undergo frequent 
collisions transporting the system collectively sidewards until the flow-induced 
pressure pushes them into the vacuum. Hence, the suppression of particle emis- 
sion at r t — > is basically caused by the dynamical expansion: many nucleons 
are transported through the medium before they reach the 'surface' and freeze 
out. Table 2 gives the average values for freeze-out radii, collective flow veloc- 
ities and transverse momenta of particles in the central rapidity region. Note 
that the average transverse 'velocities' < p t > /A, the transverse flow veloci- 
ties and source radii decrease with increasing A and saturate for cluster masses 
A > 3. The relative suppression of cluster formation at the 'surface' (r t > 6fm) 
contradicts the present fireball analyses which assume a common density and 
velocity profile for all particles. 

Clusters are clearly dominated by the collective flow components in the fi- 
nal phase space distribution: The collective transverse velocities of heavier mass 
clusters already provide most of the total transverse momentum («80%). There- 
fore, the freeze-out density and collective velocity profiles determine the final 
spectra almost exclusively: 

It is convenient to approximate the collective velocity profile by f3t—A{r t — 
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ro) B in order to demonstrate the effect of the interwining of collective velocities 
and freeze-out probabilities on the transverse spectra. This parametrization 
leads to a transverse momentum spectrum of the form 

dN _ m\ K_ TO 2)(i-B)/s 



p t d Pt BA 2 / B K)(i+b)/b 



p(r t ) ; m t :=Jm 2 a +p 2 t (22) 



where p—1/ (rt — ro)dN/dr t . The only quantities which determine such a purely 
flow-induced spectrum are the shape of the velocity profile defined by the pa- 
rameter B and the density distribution p. Assuming a box density profile the 
spectra have a convex shape and maximum at finite p t only for B<1 - consistent 
with the RQMD results. A quadratic r t -dependence (B—2) would instead yield 
an overall concave spectrum, in particular diverging at pt — > 0. 

The flattening of the transverse spectra at low pt-values is due to the suppres- 
sion of clusters in the very central region of the reaction (r t — > 0). The different 
apparent temperatures at high p t -values are caused by the strong weight of 
large flow velocities for r t > 6fm. The peak/shoulder in the transverse spectra, 
however, appears approximately at pt/A w< (it > and directly measures the 
strength of the transverse flow at the position where most of the clusters at cen- 
tral rapidities freeze out (r t ~ 5 — 7fm). Note that it is not possible to describe 
the transverse spectra with one single temperature and collective flow velocity in 
contrast to what has been claimed for reactions in the lGcV/n regime |3?], ItSSf . 

Si+Si collisions do not provide a transverse expansion comparable to massive 
reactions (Fig. |[ Table 2). Most of the clusters are emitted close to the beam 
axis where the transverse flow velocities are small. In fact, the 'surface suppres- 
sion' acts more strongly in the case of the smaller system, in accordance with 
the larger surface to volume ratio. Note that here the average transverse flow 
velocity of 4 He states is approximately a factor of two smaller than for protons 
(Table 2). Transverse flow is nearly invisible due to large 'local' momentum 
components. This sampling of clusters at smaller distances from the beam axis 
explains why the transverse flow features are mostly invisible, although the ve- 
locity profiles for cluster states in Fig. || - the 'collective' position-momentum 
correlations - are almost equal for Si+Si and Au+Au collisions. 

The role of the shapes of collective velocity and density profiles has been the 
subject of much previous debate [|l7[ [2^, [H ^38|. In particular the low-pt 
pion enhancement and the spectra of protons and deuterons have been inter- 
preted in terms of transverse flow with the assumption of an expanding thermal 
fireball [p9l [701]. A similar picture has been used to extract 'local temperatures' 



and chemical potentials exploiting the final particle ratios |71j, 72 . We wish to 
discuss this issue in more detail because in those analyses the freeze-out profiles 
used differ significantly from the results of the transport calculations presented 
here. 

In Ref. p9[ | a value for B = 2 in combination with a box-shaped density 
profile was used to explain the low-pt pion enhancement i.e. a concavely curved 
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pt-spectrum. As a consequence of these assumptions proton and deuteron spec- 
tra show the same behaviour, in particular for pt — > 0. The microscopic calcula- 
tion, however, contradicts this picture and shows profiles which are compatible 
with B Ri 0.5 and a non-trivial position geometry. The concave curvature of 
transverse spectra in the transport calculations has been confirmed by recent 
data (see protons and deuterons in Si+A reactions at 14.6AGeV and the 
results for Au(11.6AGeV)Au in Fig. |J). In this sense, neither the assumption 
p = const nor a shape of the freeze-out profile according to B = 2 -as used in 
[[39)- can be justified. The main reason for the misleading results in [|6S|] is prob- 
ably the misinterpretation of concavely shaped pion spectra. Pions are strongly 
influenced by the final decays of resonances such as A, p, B* (see Q and Refs. 
therein). The alternative prediction that the low-pt pion excess at AGS energy 
comes from A-resonances |74| fto| has been confirmed by experimental recon- 
struction of the p7r invariant masses which show a strong A-signal, in agreement 
with RQMD Furthermore, the early preliminary data for protons used in 
[[39] were limited in acceptance (m t — mo > 200MeV). They excluded those 
regions where most of the shoulder-arm effect appears and were - within the 
error bars - also consistent with concavely shaped distributions for protons and 
deuterons. 



5.2 Directed Flow 

Besides the characteristic signals in the inclusive spectra, the correlation be- 
tween rapidity and directed transverse momentum p x (y) in Fig. |] is another 
indicator of a non-trivial event geometry. This observable is well known as 
the nuclear "bounce-off" discovered first at the Bevalac J7f|. The averaged 
transverse velocity p x /A in Fig. ^| is defined by the averaged transverse momen- 
tum per cluster nucleon projected on the theoretical reaction plane for particles 
within a certain rapidity interval Ay: 



p x (yo)/A := / — ^ e x ■ fi/A 



(23) 

\y-yo\<Ay/2 



e x is the unity vector which points perpendicular to the beam axis into the 
impact parameter plane. (...) denotes the final event averaging. In the case 
of strong longitudinal and transverse flow contributions in the final source this 
quantity reflects the collective sideward flow of matter predicted by hydrody- 
namics [n~5L [76| and microscopic models JTT 
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Clusters exhibit larger p x {y)/A values than nucleons although the division 
by A excludes the trivial effect of the momentum scaling with mass pa/ A ~ p^ 
at equal velocity. This stronger correlation for cluster states is well known from 
Au+Au reactions in the lGeV/n energy regime [jf?], EHI. The reason for this 
behaviour is demonstrated in Fig. |?] which shows profiles for the collective in- 
plane velocity (3 X = e x ■ p/E and freeze-out density of protons and deuterons. 



13 



Only particles in the forward hemisphere 2.1 < y < 2.6 are taken into account. 
The profiles are drawn as a function of the transverse distance to the beam axis 
taken in the original nucleon-nucleon cms and projected onto the theoretical 
reaction plane x := x ■ e x . 

The densities for deuterons are scaled by an arbitrary factor to exhibit the 
qualitative difference between nucleons and deuterons: In contrast to the aver- 
age values, the 'local' velocities of protons and deuterons are equal. The density 
distribution of the deuterons, however, exhibits a shift towards the outward re- 
gions as compared to protons. This suppression of cluster formation near to the 
original beam axis is caused by higher relative momenta for nucleons. Therefore, 
the high transverse in-plane velocities are more strongly weighted in the case of 
cluster formation which leads to higher average velocities. Note the qualitative 
difference between this increase in the reaction plane in contrast to decreasing 
values for < pt > /A at central rapidities. 

5.3 Mean Field Dependences 

The results in calculations for Au+Au show higher longitudinal (Fig. ||) and 
transverse momenta (Fig. Table 2) caused by the additional pressure which 
is built up by the repulsive mean fields at high baryon density (up to 8po is 
achieved @). Note that the region of highest compression (p/p >3) is large 
(Y ^several hundred fm 3 ) and contains up to 60% baryons in resonance states 

The difference between potential and cascade calculations is largest in the 
low-p t part of the spectra (Fig. For heavier clusters the distributions close 
to pt = change by up to a factor of three. Nuclear matter at midrapidities 
is mostly affected by the mean field contributions at high baryon densities: 
Figs. H and || show the rapidity dependence of the proton and deuteron yields 
at low transverse momenta - here defined by a cut in transverse momentum 
Pt/A < 0.5GeV - and the average transverse momenta of p, d and 4 He. While 
cascade calculations exhibit a clear peak in the dN/ dy spectra, the calculations 
including potentials show a dip even for central events. This is due to the 
additional longitudinal expansion caused by the baryon potential interaction. 
The potentials also change the average transverse momenta by <20%. 

For the smaller system Si+Si the average transverse momenta change by less 
than 3% including potential-interaction. In earlier work we have shown that 
the life time of the reaction is not long enough to establish a thermalized high 
density phase |l4| . The 'transverse communication' is in light systems much 
smaller and does not allow for a considerable transverse push due to the mean 
fields, although, potentials also lead to qualitative changes in the distribution 
of the longitudinal momenta for the small system: The results in Fig. || show 
that cascade calculations exhibit a concave shaped spectrum which turns to 
a convex distribution if potentials are included. As in Au+Au reactions the 
yields of cluster states are most sensitive: The dN/dy values for 4 He states at 
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midrapidity differ by almost a factor of 2 between the cascade and the potential 
calculations. 

The value of the flow-correlation p x (y) in central Au(11.6AGeV)Au collisions 
is roughly a factor 1.5-2 higher due to the additional sideward push of the mean 
fields (Fig. ^J). In recent work we have shown that pions and other produced 
hadrons (antikaons, antinucleons) show a characteristic 'anti-flow' |Q caused 
by scattering off spectator-like matter. This behaviour has previously been 
discussed in reactions at lAGeV |o[ |6^| . Fig. || includes this pionic anti- 
flow which also appears to be sensitive to the baryonic mean fields, too: While 
cascade calculations show sizeable p x /m values for 7rs, the inclusion of baryonic 
potentials leads to almost vanishing p x (y)-values in the laboratory frame. In the 
principal axis system, however, the strong anticorrelation of pions to baryons 
is conserved. In the work of Li et al. |8lJ the in-plane pion flow has been 
investigated in the framework of the cascade model ART. In these calculations 
the sign of the pionic flow is equal to the baryon flow in central events which 
is in qualitative difference to the RQMD results. Note the strong dependence 
of the in-plane pion flow on different absorption rates at high baryon densities 
which has recently been analyzed for reactions at lower incident beam energies 
with the QMD-model @. 

Both results for central Au(11.6AGeV)Au collisions including baryon poten- 
tials (high in-plane flow for nucleons, small p x /m for pions) are in quantitative 
agreement with preliminary flow measurements from E877 ]82] | and E866. The 
convex proton rapidity distribution in Si+Si and the width of the distribution 
in Au+Au reactions, including potential interactions, are also in accord with 
published |6l[] and preliminary [2^] data. Nevertheless, even potential cal- 
culations overestimate stopping and underestimate the transverse momentum 
production at forward rapidities (y > j/atjv) in asymmetric reactions like Si+Au 
[psj. The results for deuterons in Fig. ^ show the same trend. 

6 Summary and Outlook 

We have presented a phase space coalescence model for cluster formation with 
A > 2 in relativistic nucleus- nucleus collisions using the Wigner-function method. 
The combination with the transport approach RQMD allows the calculation of 
cluster momentum distributions. Specifically, we have studied central Au(11.6- 
AGeV)+Au and Si(14.6AGeV)+A reactions. The microscopic model shows that 
strong stopping results in observable collective behaviour of the stopped baryon- 
rich matter. Considerable flow (< >«0.5 c) develops due to the internal 
pressure of the dense matter. The transverse expansion is most visible in the 
momentum spectra of nuclear clusters which deviate markedly from thermal 
distributions. In central Au(11.6AGeV)Au collisions the transverse spectra ex- 
hibit a strong shoulder-arm shape which is most prominent for heavier mass 
clusters. 
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The 'apparent' temperatures at high transverse momenta, which result from 
the overlay of rather small local momentum fluctuations and collective flow ve- 
locities, increase with cluster mass. Furthermore, a clear "bounce-off" event 
shape is seen in massive reactions like Au(11.6AGeV)Au. The averaged trans- 
verse flow velocities in the reaction plane < p x > (y)/A are markedly larger 
for clusters than for protons (by a factor of 2). Both the shoulder-arm shape 
and the large bounce-off signal for nuclear clusters are directly related to the 
freeze-out geometry and flow correlations. In contrast to the directed flow our 
results show smaller freeze-out radii and smaller average transverse 'velocities' 
< Pt > I A for clusters at midrapidities than for free protons. This result is in 
contradiction to simple fireball and blast wave models. 

The cluster spectra and the in-plane flow may change markedly if baryon 
potential-interactions are included: According to the calculations presented here 
the yields of nuclear clusters decrease by up to a factor of three, at low p t and 
central rapidities, in the reaction Au(11.6AGeV)Au. The average pt/A val- 
ues are ~15% harder in calculations with potentials than in the cascade mode. 
The p x (y)-correlation for nucleons and nuclear clusters increases by a factor 
of 1.5-2 while the anticorrelated in-plane flow of pions vanishes. For central 
Si(14.6AGeV)Si reactions the potentials play a negligible role in transverse di- 
rection, but affect the proton and cluster rapidity spectra. Cascade calculations 
exhibit concave spectra which turn into convex distributions if potentials arc 
taken into account. The absolute yield for 4 He clusters changes by almost a 
factor of two at midrapidity. The strong sensitivity of nuclear clusters to the 
collective flow encourage the quantitative study of the transient pressure in 
nucleus-nucleus collisions. 
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Deuteron Formation 
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Figure 1: Rapidity distributions of deuterons in Si+Si (b< lfm), Si+Cu (b< 
1.5fm) and Si+Au (b< 3fm) reactions at 14.6AGeV calculated from RQMD 
simulations including potential interactions for baryons (solid histograms). The 
symbols show E802-data from Ref. fnl for central Si+AI, Cu and Au reactions. 
Note that the data have been extrapolated in m* and contain w 15% systematic 
uncertainty. 
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Figure 2: Transverse mass spectra for protons and deuterons in central 
Au(11.6AGeV)Au reactions at ulab = 1-3. RQMD-simulations including po- 
tential interactions for baryons (histograms) are compared with preliminary 
E866-data (symbols) @. 
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Figure 3: Rapidity distributions for p,d,t and 4 He in Au(11.6AGeV)Au, b< 3fm 
(a) and Si(14.6AGeV)Si, b< lfm (b). Calculations with baryon potentials are 
denoted by bold solid histograms. Cascade calculations are shown by thin solid 
histograms. The inclusion of potentials at high baryon densities leads to stronger 
longitudinal expansion in both systems. 
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Figure 4: Transverse momemtum spectra for p,d,t and He in 

Au(11.6AGeV)Au, b< 3fm (a) and Si(14.6AGeV)Si, b< lfm (b) at central 
rapidities. Calculations including baryon potentials (bold solid histograms) are 
compared with cascade simulations (thin solid histograms). The smooth solid 
lines show Boltzmann parametrizations adjusted to the high momentum part of 
the spectra (see text) in calculations with potential interaction. 



25 



Au(11.6AGeV)Au 
b<3fm 



Si(14.6AGeV)Si 
b<1fm 



10 



10"' 



10"' 



10"' 



10"' 



10"' 




I I I I I I I 


I I I I I I I 

p (c) 






' — I t 





4 He \- 





iii 


i i i i i i i i i i i 
(d) 




P 




1 














t 


4 H £ 






10" 



- 10"' 



10"' 



10"' 



10"' 



0.0 5.0 10.0 15.0 20.0 0.0 5.0 10.0 15.0 20.0 
r, (fm) r, (fm) 



Figure 5: Freeze-out profiles of protons and tritons in RQMD-calculations 
including baryon potentials for Au(11.6AGeV)Au, b< 3fm (a,c) and 
Si(14.6AGeV)Si, b< lfm (b,d) at central rapidities. The upper part shows 
transverse velocity profiles for protons and tritons. The lower part shows the 
distributions of transverse freeze-out densities of p, d, t and 4 He. 
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Figure 6: p x (y)/j4-correlations for p, d, 4 He (a,b) and pions (c) in central 
Au(11.6AGeV)Au (b< 3fm) reactions. The figure shows a factor of two increase 
of the cluster flow if baryon potentials are included. The additional rotation of 
the event plane due to the potentials leads to an apparent vanishing of the pion 
flow in the laboratory system which is, however, still pronounced in the principal 
axis system of the rotated (baryon) event. 
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Figure 7: In-plane freeze-out velocity (a) and density profiles (b) of protons 
and deuterons in central Au(11.6AGeV)Au reactions (b< 3fm). Selected are 
particles in the forward (2.1 < ?/Lab < 2.6) hemisphere, x denotes the projection 
of the freeze-out position onto the theoretical reaction plane. The deuteron 
density is scaled by a factor 16.8 to demonstrate the increase in the transverse 
freeze-out distances from the beam axis between protons to deuterons. 
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Figure 8: Comparison of dN/dy-distributions of protons and deuterons includ- 
ing a transverse momentum cut pt/A < 0.5AGeV in calculations with (bold solid 
histograms) and without baryon potentials (thin solid histograms) for central 
Au(11.6GeV)Au (a) and Si(14.6AGeV)Si (b) reactions. 
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Figure 9: Average transverse momenta < p t > (y) of p, d, 4 He in cal- 
culations with (histograms) and without (lines) baryon potentials for central 
Au(11.6GeV)Au (a) and Si(14.6AGeV)Si (b) reactions. 
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Tabic 1 


cluster 


radius (fm) 


D (GeV 3 ) 


3 H 


1.7 


0.27 x 1CT 4 


3 He 


1.8 


0.34 x 10~ 4 


4 He 


1.5 


0.53 x 10~ 4 



Tab. 1 Root mean square charge radii for nuclear clusters with A=3,4 and 
the corresponding coupling strength in the harmonic oscillator approach. 
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Tabic 2 




n 


P 


d 


t 


a He 


4 Hc 


n 


P 


d 


t 


a Hc 


4 Hc 


Au+Au cascade 


Au+Au potentials 


<Pt> 


.5 


.5 


.42 


.37 


.38 


.36 


.54 


.54 


.47 


.43 


.44 


.42 


^J<rf> (fm) 
< p t > /A(GeV/c) 


9.5 


9.7 


7.9 


6.6 


6.2 


6.0 


10.0 


9.9 


7.9 


6.6 


6.9 


6.2 


0.71 


0.71 


0.51 


0.44 


0.44 


0.43 


0.77 


0.78 


0.61 


0.51 


0.53 


0.50 


Si+Si cascade 


Si+Si potentials 


<Pt> 


.39 


.39 


.32 


.26 


.26 


.24 


.41 


.40 


.31 


.24 


.26 


.23 


\J <r\ > (fm) 


3.9 


3.7 


2.7 


2.2 


2.2 


1.8 


3.9 


3.9 


2.7 


2.2 


2.3 


2.1 


< p t > /A (GcV/c) 


0.62 


0.61 


0.42 


0.32 


0.33 


0.29 


0.64 


0.63 


0.43 


0.33 


0.34 


0.30 



Tab. 2 Average transverse freeze-out velocities < Pt >, freeze-out radii 
\J <r\ > and average transverse momenta < p t > /A for nucleons and cluster 
with A < 4 in central Au(11.6AGeV)Au and Si(14.6AGeV)Si reactions. The 
table contains calculations with (r.h.s.) and without (l.h.s.) baryon potential- 
interaction. Only particles in the midrapidity region 1.4 < y < 1.8 (Au+Au) 
and 1.6 < y < 1.8 (Si+Si), respectively, are taken into account. 
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